Critical rotation of a harmonically trapped Bose gas 
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We study experimentally and theoretically a cold trapped Bose gas under critical rotation, i.e. 
with a rotation frequency close to the frequency of the radial confinement. We identify two regimes: 
the regime of explosion where the cloud expands to infinity in one direction, and the regime where 
the condensate spirals out of the trap as a rigid body. The former is realized for a dilute cloud, 
and the latter for a Bose-Einstein condensate with the interparticle interaction exceeding a critical 
value. This constitutes a novel system in which repulsive interactions help in maintaining particles 
together. 
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The rotation of a macroscopic quantum object is a 
source of spectacular and counter-intuitive phenomena. 
In superfluid liquid helium contained in a cylindrical 
bucket rotating around its axis z, one observes the nucle- 
ation of quantized vortices for a sufficiently large rotation 
frequency f2 0] . A similar phenomenon occurs in a Bose- 
Einstein condensate confined in a rotating harmonic trap 
@j @; D • I n particular, vortices are nucleated when the 
rotation resonantly excites surface modes of the conden- 
sate. This occurs for particular rotation frequencies in 
the interval < Q < ujj_ / \/2, where luj_ is the trap fre- 
quency in the xy plane perpendicular to the rotation axis 
z. 

Several theoretical studies have recently considered the 
critical rotation of the gas, i.e. O ~ oj±, which presents 
remarkable features §, |, §, §, £(J || [If, || @ . From 
a classical point of view, for il = lo±_ the centrifugal force 
compensates the harmonic trapping force in the xy plane. 
Hence the motion of a single particle of mass m in the 
frame rotating at frequency Q is simply due to the Cori- 
olis force 2mr x ft. This force is identical to the Lorentz 
force acting on a particle of charge q in the magnetic 
field B = 2(m/q) il. The analogy between the motion of 
charged particles in a magnetic field and neutral particles 
in a rotating frame also holds in quantum mechanics. In 
this respect, a quantum gas of atoms confined in a har- 
monic trap rotating at the critical frequency is analogous 
to an electron gas in a uniform magnetic field. One can 
then expect ||, [lO| to observe phenomena related to the 
Quantum Hall Effect. 

This paper presents an experimental and theoretical 
study of the dynamics of a magnetically trapped rubid- 
ium ( 87 Rb) gas stirred at a frequency close to ui±. We 
show that the single particle motion is dynamically un- 
stable for a window of frequencies O centered around 
u>±. This result entails that the center-of-mass of the 
atom cloud (without or with interatomic interactions) 
is destabilized, since its motion is decoupled from any 



other degree of freedom for a harmonic confinement. This 
also implies that a gas of non-interacting particles "ex- 
plodes" , which we indeed check experimentally. When 
one takes into account the repulsive interactions between 
particles, which play an important role in a 87 Rb con- 
densate, one would expect naively that this explosion is 
enhanced. However, we show experimentally that this 
is not the case, and repulsive interactions can "main- 
tain the atoms together". This has been predicted for 
a Bose-Einstein condensate in the strongly interacting - 
Thomas- Fermi (TF)— regime || . Here we derive the min- 
imal interaction strength which is necessary to prevent 
the explosion. This should help studies of the Quantum 
Hall related physics in the region of critical rotation. 

Consider a gas of particles confined in an axisymmet- 
ric harmonic potential Vo (r) , with frequency lo z along the 
trap axis z, and u± in the xy plane. To set this gas into 
rotation, one superimposes a rotating asymmetric poten- 
tial in the xy plane. In the reference frame rotating at 
an angular frequency Q around the z axis, this potential 
reads Vi(r) = emuj^Y 2 - X 2 )/2, where e > 0. The ro- 
tating frame coordinates X, Y are deduced from the lab 
frame coordinates x, y by a rotation at an angle fit. 

For a non-interacting gas, the equation of motion for 
each particle reads: 

X - 2QY + (uo]_{l - e) - fl 2 ) X = (1) 
Y + 2QX + (ujl(l + e) -fl 2 )Y = 0, (2) 

while the motion along z is not affected by the rotation. 
One deduces from this set of equations that the motion 
in the xy plane is dynamically unstable if the stirring 
frequency fl is in the interval [u>± yl — e, uj±_ \/l + e] . In 
particular, for 17 = uj±_ and e <C 1, one finds that the 
quantity X + Y diverges as exp (eu>±t/2), whereas X — Y 
remains finite. 

To test this prediction we use a 87 Rb cold gas in a Ioffc- 
Pritchard magnetic trap, with frequencies uj x = lu v — 
2ir x 180 Hz, and cu z = 2tt x 11.7 Hz. The initial tern- 
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FIG. 1: Center-of-mass displacement after free expansion 
(log-scale) vs. stirring time for fl — cu_i_ and e = 0.09. (a) Non- 
condensed cloud with 10 7 atoms, T — 5 ^K; (b) Condensate 
with 10 5 atoms. Solid line: exponential fit to the data. 



perature of the cloud pre-cooled using optical molasses 
is 100 /iK. The gas is further cooled by radio-frequency 
evaporation. For the first set of experiments we stop 
the evaporation before the Bose-Einstein condensation is 
reached. The resulting sample contains 10 7 atoms at a 
temperature T ~ 5 /iK. It is dilute, with a central density 
~ 10 12 cm -3 , and atomic interactions can be neglected 
(mean- field energy <C ksT). The second set of experi- 
ments corresponds to a much colder sample (T < 50 nK) , 
i.e. to a quasi-pure condensate with 10 5 atoms. 

After evaporative cooling, the atomic cloud is stirred 
during an adjustable period t by a focused laser beam 
of wavelength 852 nm and waist wq = 20/im, whose po- 
sition is controlled using acousto-optic modulators ||. 
The beam is switched on abruptly and it creates a rotat- 
ing optical-dipole potential which is nearly harmonic over 
the extension of the cloud. We measure the transverse 
density profile of the condensate after a period of free 
expansion. In this pursuit, we suddenly switch off the 
magnetic field and the stirrer, allow for a 25 ms free-fall, 
and image the absorption of a resonant laser by the ex- 
panded cloud. The imaging beam propagates along the 
z axis. We fit the density profile of the sample assum- 
ing a Gaussian shape for the non-condensed cloud, and 
a parabolic TF shape for the quasi-pure condensate. We 
extract from the fit the long and short diameters in the 
plane z = 0, and the average position of the cloud. The 
latter gives access to the velocity of the center-of-mass of 
the atom cloud before time of flight. 

The center-of-mass displacement as a function of the 
stirring time t is shown in Fig. [|. We choose here e = 0.09 
and fl = u>±, so that the motion predicted by Eqs.(|l|j|) 
is dynamically unstable. To ensure reliable initial con- 
ditions, we deliberately offset the center of the rotating 
potential by a few micrometers with respect to the atom 
cloud. We find the instability for the center-of-mass mo- 
tion both for the non-condensed cloud (Fig. |Tja) and for 
the quasi-pure condensate (Fig. |l|b). The center-of-mass 
displacement increases exponentially, with an exponent 
consistent with the measured e. 

We consider now the evolution of the size of the atom 
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FIG. 2: Long (■) and short (o) diameters of the atom cloud 
vs. stirring time, for Q — uj±. (a) Non-condensed cloud ; (b) 
Quasi-pure condensate (same parameters as in Fig. [j]). 



cloud as a function of t (Fig. |J) . The non-condensed cloud 
exhibits the behavior expected from the single particle 
dynamics, i.e. the "explosion" in the X — Y direction. 
The cloud becomes more and more elliptical in the xy 
plane. The long radius increases with time, while the 
short one remains approximately constant (Fig. ||a). On 
the opposite, we find that the condensate remains circu- 
lar (Fig. |^b), with no systematic increase in size. We then 
obtain the following counter- intuitive result: for a signifi- 
cant repulsive interaction the atoms remain in a compact 
cloud, while they fly apart if the interaction is negligible. 
We observe this stability of the shape of the condensate 
rotating at the critical velocity for e < 0.2. Above this 
value of e we find that the atomic cloud rapidly disinte- 
grates. For e » 0.3, after a stirring time of 50 ms, we 
observe several fragments in the time-of-fiight picture. 

We now perform a theoretical analysis of how the in- 
terparticle interaction stabilizes a rotating condensate. 
To this end we use the 2D (x, y) time-dependent Gross- 
Pitaevskii (GP) equation for an idealized cylindrical trap 
(bj z = 0). In the rotating frame the GP equation reads: 



l r 



-A+(l-€)X 2 +(l + e)Y 2 +2g\yj\ 2 -2nL yj, (3) 



where L is the z component of the angular momentum 
operator. In Eq.(||) the coordinates are given in units of 
the initial harmonic oscillator length \J h/mujj_, and the 
frequencies in units of luj_. The condensate wave func- 
tion tjj(X,Y,t) is normalized to unity, and the effective 
coupling constant is g = 4iraN, with a being the positive 
scattering length, and N the number of particles per unit 
axial length. The effective coupling g depends on density 
and characterizes the ratio of the mean-field interparticle 
interaction to the radial frequency lo± . 

Since the trapping potential is harmonic, the average 
center of mass motion of the condensate is described by 
the classical equations (|l||2|) and is decoupled from the 
evolution of the condensate wave function in the center 
of mass reference frame. We shall therefore restrict to 
wave functions ip centered at x = y = for all times. 

We start with a variational analysis of the steady state 
of the condensate in the rotating frame, using a Gaussian 



ansatz for the condensate wave function IfL 511 : 



ip(X, Y) <x exp(iaXY - (3X 2 /2 - jY 2 /2) 



(4) 



We use the symmetry properties of the Hamiltonian and 
assume that the condensate wave function remains invari- 
ant under the combination of a time reversal (jp — > ip*) 
and a reflection with respect to the XZ plane. This im- 
plies that the parameters a, [3, 7 are real numbers. We 
extremize the GP energy functional with respect to these 
parameters. Extremizing with respect to the phase pa- 
rameter a gives a = ^(7 — /3)/(7 + /3). As /3 and 7 should 
be finite and positive this sets the constraint a 2 < f2 2 . 
Extremizing over /? and 7 and expressing j3/"f in terms 
of a, we obtain a closed equation for a: 

(e/fi) [a 2 + 2an(n 2 - l)/e+ tt 2 } 

-{g/2ir)y / l-a 2 /n 2 [a 3 + (1 - 2£l 2 )a ~ fie] = 0. (5) 

In the non-interacting case (<? = 0) the Gaussian ansatz 
is exact, and a is the root of a quadratic equation. This 
ansatz also captures the scaling properties of the rotating 
condensate in the regime of strong interactions. In the 
TF limit (g -> 00) the first line of (§) can be neglected 
and we recover the cubic equation for a derived in |^] . 

For (? = the parameter a is complex in the interval of 
rotation frequencies \/l — e < £1 < ^/l + e, and there is 
no steady state solution for the condensate wave function 
at these fl |)| . For a finite g the lower border 0_ of this 
frequency interval remains equal to VI — e irrespective 
of the value of g. The upper border, dashed curves in 
Fig. 3, decreases with increasing g at a given e. For small 
anisotropy e < 1/5 it reaches the lower border at a critical 
coupling strength. For larger g the steady states exist 
at any f2. If e > 1/5, the upper border never reaches 
fi_ = \/l — e and for any g one has an interval of 
where steady state solutions are absent. The e = 1/5 
threshold was derived for the TF limit in |j| . 

We now analyze the time evolution of the condensate 
after the stirring potential has been switched on. For this 
purpose we use an approximate scaling approach to the 
solution of Eq.(||). We assume (and later on check) that 
the evolution of the condensate shape is well described 
by dilations with factors b u (t) and b v (t) along the axes X 
and Y , rotating at an angular frequency 0(t) with respect 
to the laboratory frame x,y. To determine b u (t),b v (t), 
and 4>(t), we write the wave function as 

1>(X, Y, t) = (b u b v )-^ 2 X (u, v, t) exp{i$(A > , Y, t)} , (6) 

where we have set u = X /b Ul v = Y/b v , and 

Y, t) = a(t)XY+ (b u /2b u )X 2 + {b v /2b v )Y 2 , (7) 

with a = -0tanh£ and = \n(b v /b u ). Then the GP 
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FIG. 3: Solid curves: upper and lower (fi-) borders 

of the instability region vs. g, for an abrupt switching of the 
rotation with e = 0.05 (a) and e = 0.25 (b). Dashed curves: 
upper border of the region where stationary states are absent. 



equation reduces to the following equation for x( u i v i ^) : 



cosh£ 



(ud v - vd u ) x = 



d 2 d 2 
2bJ,~2ti, 



9\X\ 



+ 2 (v 2 u b 2 u u 2 + v 2 v b 2 ,v 2 ) + b h 
The "frequencies" v u and u v are given by 



(8) 



v 2 u v = 1 =F ecos(2fif - 20) + a 2 T 20a + b u<v /b u<v . (9) 

In Eq.(|^) we took into account that u,v are eigenaxes 
of the condensate, which requires the absence of terms 
proportional to uv and is provided by the relation 

e sin(2m - 2(f)) - & - a(b u /b u + b v /b v ) + <f>£ = 0. (10) 

We now replace the lhs of Eq.(||) by fix, where p, follows 
from the normalization condition for x- The solution of 
the resulting equation is a function of b UtV ,4> and v u ,v ■ We 
then require that (§) is a relevant scaling transform, i.e., 
that the function x( u , v ,t) is most similar to the initial 
function ^(u,«,0). More precisely we set the averages 
(u 2 )t, (v 2 )t equal to their values at t = 0. This fixes v u 
and v v in terms of b u , b Vl <fi. The solution of Eqs.(^|) and 
( |l0| ) then gives the desired scaling parameters. 

The omitted lhs of Eq. (||) only insignificantly changes 
v u and v v . It vanishes in both the TF regime and for 
an ideal-gas condensate. For the TF limit our procedure 
gives the same results as the scaling approach of |l(| . 

For an abrupt switching of the rotating potential we 
use the initial conditions b UtV (0) = 1, b UyV (0) = 0(0) = 
0(0) = 0. We find two types of solution: (i) oscillating 
functions b UiV (t), (ii) one of the scaling parameters even- 
tually grows exponentially. Case (ii) describes an infinite 
expansion of the condensate in one direction, similarly to 
the expansion of the ideal gas under rotation. 

For a given e we obtain the upper and lower (f2_ ) 
instability borders in the Q — g space (see Fig. [|). The 
lower border is always equal to vT— e. The upper border 
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FIG. 4: Critical coupling strength g c vs. e for = 1. The 
solid and dashed curves correspond to an abrupt and adia- 
batic switching of the rotation, respectively. 



£l+(g) decreases with increasing g and for e < 0.17 it 
reaches f2_ at a critical value of the coupling strength. 
For e > 0.17 we have Cl + > f2_ at any g. 

The obtained results can be understood by comparing 
the frequency Q q of the rotating quadrupole mode of the 
condensate with the rotation frequency f2 ~ 1 of the 
perturbation Vi(r). In the absence of interaction one 
has ft q = 1, and the corresponding resonance leads to 
the condensate explosion. The interactions reduce the 
frequency of the rotating quadrupole mode (Q q = 1/V2 
in the TF limit), suppressing the resonance at CI w 1: 
the deformation of the condensate induced by Vi remains 
small, at least for e smaller than the detuning from the 
resonance. For larger e the condensate explodes. 

In the presence of interactions (l/v2 < Q q < 1), one 
could expect naively that the explosion occurs for a reso- 
nant drive with Cl ~ Q q , even for small e. This is not the 
case because of a non-linear character of the dynamics. 
As the system starts to elongate under the action of the 
resonant excitation, it becomes closer to an ideal gas, for 
which Q q = 1. The gas is then driven away from the 
resonance and its deformation stops. This explains why 
the lower instability border fi_ is independent of g. 

The scaling method can also be used to identify sta- 
tionary solutions, by setting if) = CI and constant b UiV 
in Eqs.(p|-[lO|). The results nearly coincide with the ones 
from the Gaussian ansatz. The existence of these solu- 
tions can also be explored using an adiabatic switching 
of the rotating potential. As shown in Fig. || the domain 
of instability for an abrupt switching of the rotation in- 
cludes the domain for the absence of stationary solutions. 

In Fig.|] we present the minimum coupling strength 
g c (e) required for the stability of the shape of the con- 
densate, for both abrupt and adiabatic switching of the 
potential rotating at ft = 1. 

To summarize, our analysis shows that the condensate 
can preserve its shape and size for rotation frequencies 
in the instability window for the center of mass motion, 
VI - e <Cl< VI + e 



cussed in 
state 




tation is switched on. For TF condensates this is the 
case if e < 0.28, which explains why the repulsive inter- 
particle interaction maintains particles together in our 
experiment. On the other hand, for larger e there are 
always frequencies at which even TF condensates are un- 
stable. This gives an account for the destruction of the 
condensate in our experiment at e « 0.3. One can think 
of observing related effects in rotating ion clouds in elec- 
tromagnetic traps (see jl7| and refs. therein). 

Even though our picture properly describes the exper- 
iment, we are likely not to deal with the ground state 
of the system. On a very long time scale, the rotating 
gas can evolve to a more complex state, for instance to a 
multi-vortex state (with possible quantum melting) dis- 
" | [lj, II, 0, or to a Quantum-Hall-like 
However these states have been discussed 
for the axially symmetric case and a natural development 
will be to include a finite rotating anisotropy e, which 
is a necessary ingredient for most present experiments. 
In particular, for CI = ujj_\/1 — e, one reaches a 1-body 
Hamiltonian corresponding to an unbound motion (with 
a gauge field) in the X direction, similar to a Quantum 
Hall fluid in a narrow channel. We believe that the study 
of many body aspects of this regime will bring in new fea- 
tures of quantum mesoscopic physics. 
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